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Abstract 

The short review of theoretical aspects of ultra high energy (UHE) neutrinos. The accelerator sources, such as 
Supernovae remnants, Gamma Ray Bursts, AGN etc are discussed. The top-down sources include Topological Defects 
(TDs), Superheavy Dark Matter (SHDM) and Mirror Matter. The diffuse fluxes are considered accordingly as that of 
cosmogenic and top-down neutrinos. Much attention is given to the cascade upper limit to the diffuse neutrino fluxes 
in the light of Fermi-LAT data on diffuse high energy gamma radiation. This is most general and rigorous upper limit, 
valid for both cosmogenic and top-down models. At present upper limits from many detectors are close to the cascade 
upper limit, and 5 yr IceCube upper limit will be well below it. 
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1. Introduction 

Many neutrino telescopes, using different technique 
of observation, were searching for High Energy (HE) 
cosmic neutrinos during last 30-40 years. Among ex- 
isting neutrino telescopes there are deep underwater/ice 
detectors (Baikal, ANTARES, IceCube, NESTOR), Ex- 
tensive Air Shower detectors (Auger and HiRes) and 
radio-telescopes: ANITA-lite, RICE, GLUE, FORTE 
and others. Even bigger projects include IEM-EUSO 
and Super-EUSO. 

Why is there such great interest to HE neutrinos? 

HE neutrinos can provide us with most important in- 
formation in physics and astrophysics. 

Detection of HE neutrinos from SN remnants will 
prove that these objects are sources of galactic cosmic 
rays (CR) and the Standard Model of Galactic Cosmic 
Ray origin will be confirmed. 

Jet models of Gamma Ray Bursts (GRBs) and Active 
Galactic Nuclei (AGN) can be proved. 

Detection of cosmogenic neutrinos can clarify the 
origin of Ultra High Energy Cosmic Rays (UHECR) 
and determine the model of transition from galactic to 
extragalactic CRs. 



Detection of neutrino jets with energies above 10 - 
10 21 eV means discovery of Topological Defects, im- 
portant objects in standard cosmology. 

Registration of HE neutrinos from the center of the 
Sun or Earth indicate the annihilation of Dark Matter 
(DM) particles there. 

Mirror matter can be discovered with help of oscilla- 
tion mirror neutrinos into visible ones. 

2. Basics of HE neutrino astrophysics 

We summarize here some basic features of HE neu- 
trino astrophysics. 

Production of UHE cosmic neutrinos 
occurs in pp and py collisions of UHE protons with the 
target nuclei and with low-energy photons. They can be 
also produced by annihilation of DM particles and by 
decays of superheavy particles. In all these cases neu- 
trinos are produced in the chain of pion decays. 

UHE neutrino sources 

are subdivided into accelerator and top-down sources, 
where neutrinos are produced in decays and annihilation 
of heavy particles. The examples of such sources are 
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given by annihilation of neutralinos in the Sun and Earth 
1 1 1, by topological defects (for a review see 0), which 
produce superheavy unstable particles, and by decays of 
quasi-stable superheavy DM particles |fj)- 

Usually HE neutrinos are accompanied by other ra- 
diations, most notably by HE gamma-rays and cos- 
mic rays (CR). There are, however, so called "hid- 
den sources" where all accompanying radiations are 
strongly or fully absorbed. The examples of such ob- 
jects are the Sun and Earth, in center of which neutrali- 
nos annihilate. Another ideal example is given by mir- 
ror matter, where all mirror particles interact with visi- 
ble matter gravitationally, and only mirror neutrinos can 
oscillate into visible ones. The almost "hidden" source 
is given by the Stecker model [4 1 of AGN, where UHE 
photons and protons are mostly absorbed or confined, 
and only HE neutrinos emerge from there. 

In Q a hidden neutrino source is produced in evo- 
lution of stellar cluster at its contraction. Collisions of 
neutron stars at the center of a cluster produce a rarefied 
cavity filled by ultra-relativistic fireballs from colliding 
neutron stars. The cavity is surrounded by thick gas en- 
velope produced by destructed stars. All HE particles, 
such as protons, gamma and X-rays are absorbed in the 
thick envelope, and only UHE neutrinos escape. 

Recently, the ultra-relativistic jet surrounded by the 
envelope of collapsing star is considered [|6 | as hidden 
UHE neutrino source. 

There are also some other more conventional exam- 
ples of astrophysical hidden sources JTHSJ. 

Neutrino detection 

includes four remarkable reactions: 

Muon production + N — > p + all gives an excellent 
tool to search for the discrete sources, since directions 
of UHE muon and neutrino coincide. 

Resonant production of W-boson, v e +e — > W~ — > v^+p 
(the Glashow resonance [9|) and v e + e — > W~ — > 
hadrons IfTUl have the large cross-sections. For the prac- 
tical applications the latter reaction is more important 
since it results in production of monoenergetic showers 
with energy E = m 2 w /2m e = 6.3 x 10 6 GeV [TO) and 
can be observable in IceCube and future 1 km 3 under- 
water detectors. 

Tau production in a detector, v T + N — > t + hadrons, is 
characterized by time sequence of three signals ifTTl : a 
shower from prompt hadrons, the Cherenkov light from 
t and hadron shower from r-decay. UHE v T are ab- 
sorbed less in the Earth due to regeneration: absorbed 
v T is converted into r, which decays producing v T again. 

Another remarkable phenomenon produced by UHE t- 



neutrinos is Earth-skimming effect [12], due to which 
the Auger observatory obtained the upper limit on UHE 
t- neutrino flux [ 1 3 1 . 

Z-bursts provide a signal from the space, caused by 
the resonant Z° production by UHE neutrino on DM 
neutrino fffl . v + vdm — > Z° — > hadrons. The en- 
ergy of the detected neutrino must be tremendous: Eq = 
m 2 z /2m v ~ 10 24 eV. 

Neutrino oscillations 

play the essential role. The neutrino flavors v e and v T are 
inefficiently produced in the accelerator sources. The 
flavor oscillation <-> v e and <-> v T can equalize the 
fluxes of these neutrinos. The oscillation length L(E) is 
given by 



L{E) = 




where Am 2 are 2.4 x 1CT 3 eV 2 and 7.7 x 10~ 5 eV 2 
for atmospheric and solar neutrino oscillations, respec- 
tively. Thus oscillation length is very short for cosmic 
distances and neutrino oscillations are very efficient. If 
neutrino flux is produced by decays of pions and muons 
with ratio k + jif~ = 1, the initial neutrino flavor ratio 
is v e : : v T = 1 : 2 : 0, and the observed flavor 
ratio (after oscillation) is v e : : v T = 1 : 1 : 1 
(equipartition). This is the case of neutrinos produced 
in pp-collisions. In case of py pion production v e can 
be strongly suppressed, but flavor equipartition after os- 
cillation is approximately holds. Many cases of flavor 
ratios with different conditions (e.g. neutron decay and 
survival, muon survival etc are studied in fl5l . 

The matter neutrino oscillations in the sources and in 
the Earth can also occur (see e.g. |6)). 

HE neutrinos from early universe. 
One might think (and many did think) that large neu- 
trino fluxes can be produced at cosmological epochs 
with large red shift z, e.g. due to decay of superheavy 
particles and production by topological defects. In fact, 
this possibility is disfavored lfl6l [TTIl by absorption of 
HE neutrinos and by nucleosynthesis bound on their 
fluxes. Neutrinos are absorbed in vv collisions with big- 
bang neutrinos and horizon of observation for neutrinos 
with energy E v q (at present) is given by redshift 

Zabs = 7.9 x 10 4 (E M /1 TeV)- 1/3 . 

Neutrino fluxes produced at large z (e.g. by topological 
defects) are strongly restricted by production of D and 
3 He at the epochs after Big Bang nucleosynthesis. Neu- 
trinos cause e-m cascades and MeV photons from these 
cascades produce D and 3 He in collisions with 4 He nu- 
clei. 
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3. Astrophysical (accelerator) sources 

Protons are assumed to be accelerated mostly by the 
shocks and produce neutrinos in pp and py collisions. 

HE neutrinos from SN remnants 

Detection of HE neutrinos from SN remnants is one of 
the most important tasks of HE neutrino astronomy, and 
this task looks perfectly realistic for IceCube and other 
1 km 3 detectors. Acceleration of protons and nuclei in 
SN remnants (SNRs) is the basic element of the Stan- 
dard Model for Galactic CR (see e.g. [18)). These 
sources successfully explain the spectra and fluxes of 
CR observed in our Galaxy, and the knee is interpreted 
as the end of Galactic CR. Accelerated protons interact- 
ing with the gas in a SNR must emit gamma-rays and 
neutrinos through production and decays of neutral and 
charged pions. HE gamma rays are detected from sev- 
eral SNRs, but from most of them the signal is compat- 
ible with the bremsstrahlung or Inverse Compton pro- 
duction by HE electron. There are in some cases the 
indications to the hadronic gamma-rays (i.e. produced 
by neutral pions born in /^-collisions), however in all 
these cases one still may argue in HE electron produc- 
tion of the observed signal. Detection of HE neutrino 
signal gives unambiguous prove of pion production by 
accelerated protons. 

There are at least two SN remnants, W28 and W44, 
for which the evidence for hadronic gamma-ray signal 
is rather strong. For both of them gamma-radiation is 
produced in a nearby dense molecular cloud. Gamma 
radiation from W28 was observed by AGILE detector 
at E > 400 MeV [19], and by Fermi-LAT J2TJ and 
H.E.S.S GUI at higher energies. The proof is based on 
spectral characteristics and molecular cloud as the tar- 
get. The case of W44 observed by Fermi-LAT is simi- 
lar: gamma-radiation is observed from dense cloud ring 
around W44 (22). 

SNRs discussed above are candidates for detectable 
neutrino fluxes by IceCube. Another candidate for the 
detected hadronic gamma-rays SNR RX J 171 3.7-3946 
was analyzed in [23 1 for detectable neutrino flux. The 
authors accurately recalculated the observed gamma- 
ray flux, assuming its hadronic origin, to HE neutrino 
flux, and found that it corresponds to 2 - 3 events per 
year at energy above 1 TeV for a detector as IceCube. 
Gamma Ray Bursts. 

GRBs are most exiting sources of UHE neutrinos. There 
are two mechanisms of HE neutrino generation. In the 
first one [24 1 particles are accelerated by external shock 
and neutrinos are produced in py collisions with GRB 
photons behind the shock. In the second mechanism 
11251 protons are accelerated by internal shocks, with 



the spectrum assumed to be oc l/E 2 . Neutrinos are 
produced in py collisions with GRB photons. All es- 
timates are very transparent and follow from assump- 
tion that the energy outputs in GRB photons, acceler- 
ated protons and produced neutrinos are about the same: 
W Y ~ W p ~ Wgrb- Then the total number of neutrinos 
with energy E per burst is 



N V (E) 



W C RB 



7-2 



In E max / E n 



where E m i n and E m3X are minimum and maximum ac- 
celeration energy, respectively. Now one can express 
the flux of neutrinos from a single GRB, in terms of 
neutrino fluence S v : 



F v + y (£) = - — 



7-2 



and calculate the number of muons produced in 1 km 3 
detector per burst as 



■.-«./ 



dE F y +y (E) o- vN (E), 



where N„ — 6 x 10 38 is the number of nucleons in 
the detector and cr vA , is v^N - cross-section. For S v ~ 
10~ 5 erg/cm 2 , with frequency of bursts as in Fermi GRB 
Monitor Nf, ~ 500 yr , one finds the number of pro- 
duced muons with E > 10 TeV for 5 yr as P/jNht ~ 0.06, 
i.e. too low for detection by IceCube. 

The diffuse flux is estimated in identical way through 
the local neutrino emissivity £. v (0) and evolutionary 
factor £ evo i 



3 An E 2 lnE m;ix /E mia 
where evolutionary factor is given by 



^evol 



dz 



/evol(z) 



(i +z) 2 V(i +z) 3 n m + A 



For the evolutionary function f evo i(z) we take the case of 
strong star-formation evolution from [26 1, which results 
in &evoi = 7.0 as a maximum value. For the emissivity 
we use 

£ y (0) < £grb(0) = 0.6 x 10 43 erg/Mpc 3 yr, 

where the local GRB emissivity is taken from 1271 . The 
most recent estimate 0.5 x 10 43 erg/Mpc 3 yr [28 1 agrees 
with the value above. For more detailed discussion see 



As a result E 2 Jy ii+Vii (E) = 1.15 eV/cm 2 s sr, i.e. 25 
times lower than sensitivity of IceCube shown in Fig|4] 
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GRB neutrinos are detectable by IceCube in case the 
hadronic energy output is an order of magnitude higher 
than one observed in photons, but it further strength- 
ens the energetic deficit in GRB models. From observa- 
tional point of view the signature of neutrino observa- 
tion is very reliable due to time and direction correlation 
with gamma-radiation. 

HE neutrinos from AGN jets. 

The models of HE neutrino production in the AGN jets 
are very similar to that for GRBs. The protons are ac- 
celerated by the multiple shocks in the AGN jets, es- 
pecially in their inner parts. Neutrinos are produced in 
the collisions with photons from the accretion disc and 
from photons produced in the jet by accelerated elec- 
trons and protons. Neutrino flux from individual AGN 
is very small, but the diffuse flux is predicted to be de- 
tectable by IceCube detector. The estimates can be per- 
formed similar to those in subsection above. One can 
found the detailed calculations in [30 3Tl l32ll . 



HE neutrinos from galaxy clusters 1331 . 

The clusters of galaxies are able to confine the UHE 
particles for a time exceeding the age of the universe. 
This is the key phenomenon which makes galaxy clus- 
ters the powerful sources of UHE neutrinos. The parti- 
cles are accelerated in clusters by various mechanisms: 
in the normal galaxies by SN shocks, in AGN and cD- 
galaxies, in the process of galactic merging etc. The 
diffuse HE neutrino flux is determined entirely by basic 
parameters characterizing the clusters. In particular, for 
the lower limit of the diffuse flux provided by normal 
galaxies in a cluster with CR luminosity L p and gen- 
eration index y g , both taken as ones in our galaxy, the 
diffuse flux is given as 

-v N g 

J V (E) oc L p E 7s — gklbPa, 

where R c i ~ 2 Mpc is the virial radius of a cluster, 
N g ~ 100 is richness of a cluster, p a is critical cos- 
mological density, and %Q.b is cosmological baryonic 
density provided by clusters. The flux is marginally 
detectable by IceCube. 



4. Non-accelerator neutrino sources 

These sources include objects with annihilation of 
DM (the Sun, Earth, cores of the galaxies), objects 
with the decays of SHDM particles (galactic halos) and 
TDs. In the last two cases neutrinos are produced in the 
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Figure 1: Diffuse all-flavor neutrino spectrum from necklaces for 
mx = 1 X 10 14 GeV (3). The thick curve gives p + y flux normal- 
ized to the AGASA UHECR data. If to normalize the proton flux in 
this figure by HiRes data all curves, including one labeled v (neutrino 
flux), should be lowered by factor 3-5. 



decays of superheavy particles with the masses up to 
Mgut ~ 10 16 GeV. Neutrino spectrum can be approxi- 
mately described at highest energies as dE/E 2 . 

Neutralino annihilation in the Sun and Earth. 

Neutralino is the best motivated DM particle. Cross- 
ing the Sun or Earth a neutralino can loose its energy 
in collisions with nuclei and diminish its velocity below 
the escape velocity. If it happens, a neutralino becomes 
gravitationally trapped in the object, and loosing further 
their energies, neutralinos are accumulated in the cen- 
ter of a celestial body [1 1. Annihilating there they pro- 
duce short-lived hadrons, e.g. D-mesons, which decay 
to neutrinos. The process of annihilation strongly de- 
pends on neutralino mass and composition (mixture of 
basic fields: zino, bino and two higgsinos). 

Superheavy Dark Matter (SHDM) 

The first proposal of SHDM [34| was motivated by Ul- 
tra High Energy Cosmic Rays (UHECR) and by natu- 
ral character of DM production at the epochs soon af- 
ter inflation. In particular SHDM particles can be pro- 
duced gravitationally 11351 . when the Hubble parame- 
ter H(t) exceeds the particle mass H( f) > mx- The ob- 
served density of DM in the universe £2 c d m ~ 0.23, de- 
termines the mass of the particle as mx ~ 10 13 GeV. 
The SHDM particles (X-particles) can be stable (due to 
e.g. discrete gauge symmetry) or quasi-stable (due to 
superweak discrete gauge symmetry breaking). The en- 
ergy spectrum of produced particles has approximately 
power-law form at the highest energies oc E~ 9 0. The 
dominant decay particles are photons and neutrinos. As 
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Figure 2: UHE neutrino flux in the dip model with AGN as the sources 
of UHECR. The cosmological evolution of AGN with m = 2.1 up to 
z c = 1.2 is taken from X-ray observations of AGN. At larger z the 
evolution is frozen up to z max = 2.0. The fit of the dip is very good, 
though requires y g = 2.52 different from the non-evolutionary case 
m = 0. The neutrino fluxes are given for one neutrino flavor. 



any cold DM, X-particles are accumulated in the halos 
of galaxies, in particular in our galaxy with overdensity 
2.1 x 10 5 . One can expect the detectable fluxes of UHE 
photons and neutrinos from the Galactic Center region. 

Topological defects (TDs). 

TDs are fundamental cosmological objects. They are 
produced in early universe due to symmetry break- 
ing accompanied by the phase transitions. In many 
cases TDs become unstable and decompose to con- 
stituent fields, superheavy gauge and Higgs bosons (X- 
particles), which then decay producing UHE neutrinos 
(see (2][36) for the reviews). 

Ordinary strings are produced by U( 1 ) symmetry break- 
ing. There are several mechanisms by which ordinary 
strings can emit HE neutrinos: collapse of the string 
loops, self-intersection, annihilation of cusps, produc- 
tion and annihilation of tiny loops. In most cases pro- 
duced neutrino fluxes are too low for detection. More 
promising scenario is given by the radiation of X- 
particles through the cusp, a peculiar point where ve- 
locity reaches velocity of light. This point appears on a 
loop during each period. The points of the loop in the 
cusp region have distribution in values of the Lorentz 
factors F from a maximum value in the cusp point to 
F = 1 . The particles escaping through a cusp region are 
boosted by these Lorentz factors. Recently, such case 
has been considered in ll37l . The bosonic (Higgs) con- 
densate in a string loop, emits Higgses through the cusp, 
and due to the Lorentz factor boost, these particles can 



reach the tremendous energies. 
Superconducting strings 

can be powerful sources of neutrinos. In a wide class of 
elementary particle models, strings behave like super- 
conducting wires. The charge carriers are massless in- 
side the string and superheavy outside. Moving through 
cosmic magnetic fields, such strings develop electric 
current. When the current reaches the critical value, the 
charge carriers escape from a string, turn into massive 
mode and decay. This process is strongly enhanced near 
the cusps due to the Lorentz boost. The decay prod- 
ucts, in particular neutrinos, are emitted isotropically in 
a frame of cusp segment at rest, and propagate in the 
laboratory system as a very narrow jet with the opening 
angle § ~ 1/F. This scenario was analyzed numerically 
in ll38ll with two main model features included. First, 
from all known structures of the universe, the excita- 
tion of electric current occurs most efficiently in clus- 
ters of galaxies for which the magnetic field reaches 
B ~ 1(T 6 G and filling factor / ~ 10 3 . Second, the 
symmetry breaking scale of order 10 9 - 10 12 GeV must 
be assumed for detectable neutrino fluxes. The typical 
Lorentz factor of the radiating cusp segment is calcu- 
lated to be r c ~ 10 12 , and the maximum energy of emit- 
ted particle can reach r c 77 ~ 10 22 GeV. The neutrino 
spectrum is assumed oc l/E 2 , similar to oc E~ 19 Q. 
The spectrum E 2 J V (E) =const can be very close to the 
£T 2 -cascade upper limit in Figj4] 

Necklaces (monopoles connected by string) 
are produced in the G — > H X U(l) — > H x Z2 sequence 
of symmetry breaking, with each monopole being at- 
tached to two strings, and a loop reminds a necklace 
with monopoles playing the role of beads. In the pro- 
cess of evolution the strings shrink due to gravitational 
radiation and MM pairs in the necklace inevitably an- 
nihilate. This model is most plausible and well devel- 
oped for UHE neutrino production. Diffuse neutrino 
flux from necklaces are shown in Fig. [T] according to 
calculations in Q. 

5. Cosmogenic neutrinos 

Starting from pioneering work ll39ll the fluxes of cos- 
mogenic neutrinos have been calculated in many works 
ll40l - |49l. The predicted fluxes differ very consider- 
ably, depending on the assumptions about mass com- 
position of accelerated particles, on maximum energy 
of acceleration and on cosmological evolution of the 
sources. We present here the UHE neutrino fluxes 
calculated in the dip model for the observed UHECR 
||29l |50l , assuming AGN as the sources (3T). This 
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Figure 3: Evolution factor k ev , which describes increasing of neutrino 
flux due to cosmological evolution of the sources. The evolution fac- 
tor is shown as function of z max for different indices y g and parameters 
of evolution m. The evolution factor k ev is large for large m and z max - 

model is valid for the proton-dominated composition of 
UHECR, based on observations of HiRes. The pair pro- 
duction dip is a feature of interaction of extragalactic 
UHE protons propagating through CMB. It is caused by 
energy losses of protons due to p + jcmb — > e + + e~ + p 
scattering. This feature in proton spectrum is well con- 
firmed by observational data Il29ll50l . 

To calculate neutrino flux produced by UHE protons 
it is enough to know the generation rate of UHE protons 
at each cosmological epoch. We take it as Q(E)(1 +z) m , 
where Q(E) cc E~ y s and (1 +z) m describes the cosmolog- 
ical evolution of the sources up to some maximal red- 
shift Zmax- In calculations we consider two cases: with- 
out evolution when we have only one free parameter, 
the generation index y g , and evolutionary scenario with 
three free parameters y g , m and z max . We must fit the 
observed HiRes spectrum with one calculated at z — 0. 
In non-evolutionary scenario the best fit with very good 
X 2 is given by y g = 2.7. For the fit with evolutionary 
model we assume AGN as the sources and take the AGN 
evolution from X-ray observations [52 1: (1 + z) m with 
m = 2.7 up to z c = 1.2, and frozen evolution from z c 
to Zmax = 2. The generation index is fixed as y g = 2.52 
for the best fit of HiRes data (see Fig. [2|. One may no- 
tice that the theoretical dip automatically describes the 
ankle at E m 5 x 10 18 eV. The calculated one-flavor neu- 
trino fluxes for the AGN evolutionary model are shown 
in Fig.|2]for two values of E m!a . 

The evolution is able to increase strongly the neutrino 
flux. The increase of the flux is given by the evolution 
factor k ev , which depends on m, z max and y g . This de- 
pendence is shown in Fig. [3] 



Figure 4: The experimental upper limits on UHE neutrino fluxes in 
comparison with e-m cascade upper limit in assumption of E gen- 
eration spectrum (curve E~ 2 cascade) and with predictions for cos- 
mogenic neutrinos in the dip model (curves dip-min and dip-AGN), 
for neutrinos from necklaces and from SHDM. Neutrino fluxes from 
necklaces and SHDM are normalized by AGASA data, and for nor- 
malization by HiRes data the fluxes should be diminished by fac- 
tor 3 - 5. Neutrino flux from superconducting strings is given by 
E 2 J(E) =const and it can reach the upper limit '£~ 2 cascade'. Neu- 
trino fluxes are given for three flavors. 

6. Cascade upper limit on diffuse neutrino flux 

The e-m cascade upper bound puts the rigorous up- 
per limit on UHE neutrino flux l53l [8) . This limit, in 
contrast to WB upper limit ll54ll . is valid for both ac- 
celerator and non-accelerator neutrinos. The produc- 
tion of neutrinos is accompanied by production of high 
energy photons and electrons from pion decays. Col- 
liding with low-energy target photons, a primary pho- 
ton or electron produces e-m cascade due to reactions 
7 + Ttar — * e+ + e i e + Ttar — > e ' + Y> etc - The cascade 
spectrum in its high-energy part is proportional to E~ 2 , 
which is very close to the EGRET observations in the 
range 10 MeV - 100 GeV ll55l . The observed energy 
density in this range is wegret ~ (2 - 3) x 10~ 6 eV/cm 3 . 
The cascade energy density must be <x> cas < w E gret> 
and it limits diffuse neutrino flux. The situation has 
dramatically changed with the new data of Fermi-LAT 
ll56l on the flux and spectrum of diffuse extragalac- 
tic gamma-radiation. In comparison with EGRET this 
flux is lower and spectrum is steeper (oc E~ 2A ). It re- 
sults in stronger upper limit on the cascade energy den- 
sity w cas < 5.8 x 10~ 7 eV/cm 3 ll57l . which severely 
diminishes the allowed UHE neutrino fluxes ll57l l58l . 
The maximally allowed cascade energy density «™ x ~ 
5.8 x 10~ 7 eV/cm 3 provides the upper limit on the inte- 
gral UHE neutrino flux J v (> E) (sum of all flavors). It 
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is given by chain of the following inequalities 

An f°° An f°° 

<y cas > — I EJy(E)dE > — El J v (E)dE 
c Je c J e 

where a> cas < <w™ x , and the integral in rhs of Eq. (|6]l 
gives the integral spectrum of neutrinos J(> E). Thus, 
this equation gives the upper limit on the integral neu- 
trino flux, which can be expressed in terms of the upper 
limit on differential neutrino spectrum J V (E) as 

E 2 J V (E) < -f <f . (1) 
An 

Eq. ([T]i gives the rigorous upper limit on the neutrino 
flux. It is valid for neutrinos produced by HE protons, 
by topological defects, by annihilation and decays of su- 
perheavy particles, i.e. in all cases when neutrinos are 
produced through decay of pions and kaons. It holds for 
arbitrary neutrino spectrum falling down with energy. If 
one assumes some specific shape of neutrino spectrum, 
the cascade limit becomes stronger. For E~ 2 genera- 
tion spectrum, which is used for analysis of observa- 
tional data one obtains the stronger upper limit. Given 
for three neutrino flavors it reads 

max 

E 2 J V (E) < ^ , (2) 

An ln(£ max /£ min ) 

This upper limit is shown in Fig. [4] 

The most interesting energy range in Fig. [4] corre- 
sponds to E v > 10 21 eV, where acceleration cannot 
provide protons with sufficient energy for production of 
these neutrinos. At present the region of E v > 10 21 eV, 
and E v » 10 21 eV is considered as a signature of top- 
down models, which provide these energies quite natu- 
rally. 

As one can see from Fig. |4]the observational upper 
limit for IceCube after 5 yr of observations will be be- 
low the cascade upper limit. Crossing it, this detec- 
tor will enter the physically allowed region of neutrino 
fluxes, and it can be regarded as historical event. The 
WB upper limit is not relevant for UHE neutrinos: it is 
not valid for top-down scenarios because proton produc- 
tion is strongly suppressed for top-down sources, and it 
is very uncertain for cosmogenic neutrinos, where for 
the same proton flux the fluxes of accompanying neu- 
trinos may differ by one-two orders of magnitudes (see 
Fig. [2). However, the WB upper bound remains the con- 
venient low-flux benchmark for detection of neutrino 
fluxes. 

Mirror matter and mirror neutrinos. 

Mirror neutrinos give the only example of fluxes not 
limited by the cascade upper limit. 



The concept of mirror matter, as first was suggested 
by Lee and Yang l59l . consists in existence of sector 
of matter fully symmetric with ordinary one and gen- 
erated by space-reflection transformation. Kobzarev, 
Okun and Pomeranchuk |60| added the basic assump- 
tion that these two sectors communicate only gravita- 
tionally. The gravitational interaction results in mixing 
of mirror and ordinary neutrinos and their oscillations 
J3T). In two-inflaton cosmological model |62l the mir- 
ror matter is suppressed, while mirror TDs can strongly 
dominate. Mirror TDs copiously produce mirror neutri- 
nos with extremely high energies, which oscillate into 
visible ones, while all other mirror particles, which ac- 
company production of mirror neutrinos, remain invis- 
ible for our detectors. Therefore, the upper limits on 
HE neutrinos in our world do not exist and their fluxes 
can be above the upper limit shown in Fig. [4] Neutri- 
nos from TDs typically have very high energies and one 
can see that fluxes of discussed neutrinos are now very 
severely constrained by ANITA-lite data |63|. 

7. Conclusions 

UHE neutrinos are expected to solve many problems 
in astrophysics and cosmology. Detection of HE neu- 
trinos from SNRs by IceCube is needed for confirma- 
tion of Standard Model for GCRs. UHE neutrinos from 
GRBs will clarify the nature of these most unusual and 
controversial objects. As far as HE radiation is con- 
cerned, this is true for AGN, too. 

The diffuse UHE neutrino radiation is presented by 
cosmogenic and top-down neutrinos, in particular neu- 
trinos from TDs. The fundamental problem of astro- 
physics involved in cosmogenic neutrinos is acceler- 
ation of particles. The shock acceleration at present 
knowledge of its theory cannot provide E mlix higher than 
10 21 - 10 22 eV, and thus energies of cosmogenic neutri- 
nos cannot exceed 3 x 10 2l) eV. TDs naturally produce 
neutrinos emitted from cusps with energies by many or- 
ders of magnitude higher. Detection of neutrinos with 
these energies mean discovery of new physics. 

Cascade upper limit is very general bound valid for 
both cosmogenic and top-down neutrinos. This upper 
limit became stronger with new Fermi-LAT data on ex- 
tragalactic HE diffuse gamma-radiation. From all exist- 
ing detectors only IceCube reached the sensitivity be- 
low the cascade upper limit (see Fig. |4]l and entered the 
physically allowed region for detectable neutrino fluxes. 
It can be considered as historical event for HE neutrino 
astronomy. 

The flux of cosmogenic neutrinos can be large only in 
case UHECR are proton-dominated. Even in this case 
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the flux is detectable if maximum acceleration energy 
Emax is large and sources have strong cosmological evo- 
lution (see 15711581 "). 

Cosmogenic neutrinos of highest energies are de- 
tectable by future experiment JEM-EUSO in rather ex- 
treme models with large E m . dx and strong cosmological 
evolution (see Il57ll58l ). 

The search for UHE neutrinos in any case is a search 
for a new physics, either for astrophysics (the new accel- 
eration mechanisms and cosmological evolution of the 
sources) or for topological defects, mirror topological 
defects and superheavy dark matter. 
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